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Abstract 

We consider an elastic rod model for twisted DNA in the plectonemic regime. The molecule 
is treated as an impenetrable tube with an effective, adjustable radius. The model is solved 
analytically and we derive formulas for the contact pressure, twisting moment and geometrical 
parameters of the supercoiled region. We apply our model to magnetic tweezer experiments of a 
DNA molecule subjected to a tensile force and a torque, and extract mechanical and geometrical 
quantities from the linear part of the experimental response curve. These reconstructed values 
are derived in a self-contained manner, and are found to be consistent with those available in the 
literature. 
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I. INTRODUCTION 



Mechanical properties of the DNA molecule play an important role in the biological pro- 
cesses involved in the cell, yet we only have an imprecise view of these properties. Advances 
in nanotechnologies make it possible to exert forces onto isolated DNA filaments: mechan- 
ical response of the molecule is now widely studied. Single molecule experiments provide a 
powerful way to investigate the behavior of DNA subjected to mechanical stress. In such 
experiments, the molecule is held by optical or magnetic tweezers and forces and torques are 
applied to it . The interaction between DNA and proteins is actively investigated; for 
instance, the chemical and mechanical action of an enzyme on a molecule can be inferred 
from the global deformation of the molecule 

In this paper we focus on a specific type of experiments: a double stranded DNA molecule 
is fixed by one end to a glass surface while the other end is attached to a magnetic bead; 
using a magnet, a pulling force and a torque are applied on the DNA filament 0]. Large 
ranges of pulling forces, from one tenth to one hundred piconewton, and number of turns 
can be explored in the experiments, and the molecule displays a variety of behaviors and 
conformations j^, [g], 3] • We study the response of the molecule to moderate forces, below 
10 pN, and moderate to large number of turns, equivalent to a positive supercoiling ratio of 
the order of 0.1. In experiments, the pulling force is kept constant while the bead is rotated 
gradually. Above a threshold value of the number of turns, the molecule wraps around itself 
in a helical way, giving rise to a structure comprising plectonemes. The vertical extension 
of the molecule is recorded and plotted as a function of the number of turns. Experimen- 
tal rotation-extension curves have a characteristic shape and are called hat curves 
The central, rounded part of the curve can be explained using the worm-like chain (WLC) 
model [l^ and its variants. At larger number of turns, the extension of the molecule de- 
creases linearly. This linear part is obtained when the molecule is in supercoiled configuration 
and forms plectonemes. The plectonemic structure is made of two interwound helical fila- 
ments whose geometry is characterized by the so-called superhelical angle and radius; note 
that each of these filaments is itself made of a double-stranded DNA molecule. The super- 
helical angle and the twis ting moment in the filaments are key parameters that control the 



action of topoisomerases 



ll| . RNA polymerase [121, other enzymes [l3| on DNA. The 



distance of self-approach of DNA in supercoiled regime has been the subject of a number 
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of studies Q, [isl, 16, 17|. In previous analytical and numerical work, the double stranded 



DNA molecule has been modelled as a twist-storing elastic filament. These approaches have 
been successful at reproducing the response of DNA to moderate torque Q, given by 
the central region of the experimental curves. The analysis of the linear regions of these 
curves, based on a detailed model of plectonemes, was lacking until recently: in Ref. [2^, a 
composite model based on an empirical free energy of supercoiled DNA is proposed. 

Here we present an elastic rod model for helical supercoiling of the DNA molecule, which is 
relevant to a large number of turns. Our model is self-contained and provides an mechanically 
accurate description of elastic filaments in contact. The molecule is divided in two domains: 
one where the configuration is a worm-like-chain, dominated by thermal fluctuations, and 
the other one, a superhelical region dominated by elasticity, where the molecule contacts 
itself. The plectonemic regions can be spread in various places of the molecule; as this 
does not change the mechanical response of the system, we refer to these regions as if 
they were in one chunk. We deal with self-contact by introducing an effective superhelical 
radius (distinct from the crystallographic radius of 1 nm, from the size of the Manning 
condensate and from the Debye length, although in the same range of values), which varies 
with external loads and salinity of the solution. The effective radius is defined as the radius 
of a chargeless, impenetrable and elastic tube having the same mechanical response as the 
molecule. This radius is not given as a parameter of the model and is extracted from 
experimental data. Using an energy approach, we relate geometrical variables (superhelical 
radius and angle) to applied force and torque. We also characterize the response of the 
molecule in the plectonemic regime, extend former numerical results 2l|], and show how 
geometrical and mechanical parameters can be extracted from experimental data. 



II. MODEL 



The present model investigates the equilibrium behavior of an elastic rod with bending 
rigidity Kq (the bending persistence length is A = Ko/^kT), where k is the Boltzmann 
constant and T the absolute temperature) and twisting rigidity under traction and 
torsion as shown in Fig. [TJ This is a coarse-grained model for DNA where base-pairs details 
are neglected. For instance, the anisotropic flexibility of the molecule, originating from base 
pairing and major-minor grove geometry, is smoothed out at a scale of several base pairs: a 
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highly twisted anisotropic rod can be replaced by an equivalent isotropic rod with effective 
bending rigidity [22 1. 



Geometry 



We start with a geometric description of the rod configurations relevant to the plectonemic 
regime. This defines a reduced set of configurations (Ansatz), over which we shall minimize 
the elastic strain energy associated with deformations. The rod, of length is considered 
inextensible and has circular cross- sect ion; let s denote the arclength along the rod. The 
strain energy involves, at lowest order, the geometric curvature k{s) of the centerline of 
the rod as well as the twist t{s). We emphasize that the twist r(s) is different from the 
geometrical (Frenet-Serret) torsion of the centerline as it takes into account the rotation of 
material cross sections around the centerline. It allows one to distinguish between twisted 
and untwisted configurations of the rod having the same centerline. The rod centerline is 
parameterized by r(s) and its unit tangent t =^ dr/ds can be described with spherical angles, 
as shown in Fig. [TJ a{s) is the zenith angle and ipi^) the azimuth angle with respect to the 
direction ex along the common axis of the two superhelices in the plectonemic region. 

We consider the following configurations, relevant to a large applied number of turns, n. 
The tails are assumed to be straight but twisted (thermal fluctuations will be accounted for 
by using the rescaled tail length predicted by WLC theory). The plectonemes are described 
by two identical and uniform helices where, again, each of these helices is itself a double- 
stranded DNA molecule. Both the end loop of the plectonemes and the matching region 
between the tails and the plectonemic part are neglected. 

In the tails the rod is straight and aligned with the e^, axis, and t = there. The 
geometric curvature k == |dt/ds| is zero, k,{s) = 0. 

In each filament of the plectonemes, the position vector r(s) and the tangent vector t(s) 
describe a superhelix of axis e^: 



The other filament of the plectonemes is obtained by a rotation of 180° around the axis 




s cos a 



cos a 




(1) 
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FIG. 1: Sketch of the magnetic tweezers experiment. A B-DNA molecule of total contour length 
I is fixed in s = to a glass surface while the other end in s = ^ is attached to a magnetic bead. 
A pulling force Fext and a torque Mext are applied at the upper end by using a magnet. The 
superhelical angle and radius are denoted a and R respectively. 



ex. Here x = stands for the chirality of the two helices and R and a denote the 
superhelical radius and angle, respectively. In equation ([1]), the condition dr/ds = t yields 



d^/'/ds = X sina/i?. The curvature in the plectonemes is =^ |dt/ds| = 

Noting £p the contour length spent in the plectonemes, we obtain the following expression 
for the integral of the squared curvature over the whole length of the rod: 

n , sin'^a „ . . 

«:M. = -^V (2) 

The end torque twists the filament. For a rod with circular cross-section, the twist r(s) 



at equilibrium is uniform 



23| . dr/ds = for all s. As a result, the equilibrium configuration 



of the rod is fully specified by the centerline, through the variables a, R and £p, and an 
additional scalar r describing twist. 

The twist parameter r is geometrically related to the number of turns imposed on the 
magnetic bead, n, which is equal to the link of the DNA molecule, n = Lk. In the present 
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case the link reads 2l| : 



Ti rr. 1 / 1 sin 2a 

Lk = Tw + Wr = — / Tds-x 



271 



AnR " 2 



sin 2a 
2R 



p ' 



(3) 



as we neglect the writhe in the tails. 



Energy formulation 



Using the above notations the rod is described by four variables: a the superhelical angle, 
R the superhelical radius, r the twist and ip the contour length spent in the plectonemes. 
We proceed to derive the total energy of the system as a function of these four variables. 
It is the sum of three terms, V = Vci + Vext + Mnt, where the first term is the strain elastic 
energy, the second is the potential energy associated with the external loads Fcxt and Mext, 
and the third accounts for interaction of the filaments in the plectonemes. The strain elastic 
energy for the rod of total contour length £ is : 



ds . (4) 

lo ^ Jo 

We do not take into account the reduction of the effective torsional rigidity in the tails due 
to fluctuations 19|]. The potential energy is given by: 



cxt 



-Fext(2;(£)-2;(0))-27rMextn, 



(5) 



where z{i) — z{0) = i — ip for straight tails and n = Lk. 



If the DNA-DNA interaction was clearly estab 



interaction energy Vint in the total energy V 2^. This is not the case and we model the 



ished, we would include the corresponding 



filaments in electrostatic interaction as effective chargeless hard-core tubes. The effective 
radius a of these tubes accounts for a variety of physical mechanisms and in particular for 



the presence of counter-ions. As in Refs. 15|, |25|, we do not fix the quantity a in advance 



and let it vary with experimental conditions, such as applied load and salinity. In fact, we 
show how a can be extracted from experimental measurements. Doing so, we replace the 
actual (unknown) interaction potential Vint{R, «) by a hard-core interaction with adjustable 
radius a, and optimize a to best fit a given experiment. 



The parameter a must certainly be larger than the crysta' 



It is different from the radius of the Manning condensate 



lographic DNA radius 1 nm. 



271, 



28| since approximately 



a quarter of the charge remains outside of the Manning condensate. The equihbrium is 
the solution of a constrained minimization problem for the elastic energy, subjected to the 
impenetrability condition 

R>a. (6) 

We anticipate on the fact that there is contact, R = a, for typical experimental conditions. 
This constraint is handled by a Lagrange multiplier A; the actual interaction energy is then 
substituted with the following expression: 

V,^t = -X{R-a). (7) 

Combining Eqs. ([2H7]), we write the total potential energy of the system as: 

Via, i?, r, Q = ^ ^ i^ + ^r'i- F,^, (i - Q 

-M^^Jri-x^^i^-XiR-a). (8) 



In Ref. 



29|] a similar energy function has been introduced but the rest of analysis differs 
from ours. Indeed, their approach focuses on statistical mechanics and the analysis of the 
state of lowest energy is overlooked. Moreover, the parameter a is fixed a priori to the 
crystallographic radius of DNA, a = 1 nm, which is a strong underestimation of the actual 
distance of self-approach of DNA in saline solution. In contrast, we undertake a detailed 
analysis of the equilibrium solutions, with thermal fluctuations considered in the tails; this 
allows us to derive simple formulas for the force and the moment as a function of the 
superhelical variables, applicable to magnetic tweezers experiments. 



III. RESULTS 



Mechanical equilibrium is given by the Euler-Lagrange condition for the stationarity of 
the potential V{a, R,T,ip) in Eq. ([8]) with respect to its variables, 

/dV_ dV_ dV_ dV\ 
\d^' 9^' di^' dRJ ~ • 

The first condition dV/dr allows one to recover the constitutive relation for twist defor- 
mations, Mext = K3 T, given that the twisting moment is uniform in the filament and equal 
to the applied torque Mext- 
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Variation of the total energy with respect to a gives the expression of the apphed torque 
Mext in terms of the superhehcal variables a and R: 

^ _2xKocosasin^a 

R cos2a ' ^ ^ 



which is what was found for purely plectonemic solution (no tails) SOj. 

The condition dV/dip = 0, combined with Eq. ([9]), allows one to relate the pulling force 
Fext to the superhehcal geometry: 

F„ = §si„*.(i + ^). JIO) 

This formula justifies and extends the numerical fit F^^t K^a^/ R^ found in Ref. |2]J for 
small values of a. 

The Euler-Lagrange condition with respect to R yields an equation involving the Lagrange 
multiplier A. The quantity A/^p can be interpreted as the contact force per unit length, p, 
of one filament onto the other. Eqs. ([8H9]), together with the condition dV/dR = 0, yields: 

^~ ip~ R^cos2a ^ ^ 

Note that this pressure (more accurately, force per unit length) is positive for a < n/A; if our 
assumption of contact R = a was incorrect, this would be indicated by a negative pressure 
value here. 

In magnetic tweezers experiments, the pulling force Fext is imposed although the applied 
torque Mext is unknown. The two unknowns R and a are then related by Eq. ( fTOj) : in the 
next Section, a second equation relating those unknowns and the extension z is given, which 
makes it possible to solve for R and a. The twisting moment can then be found from Eq. (Q. 

Vertical extension of the filament 

In magnetic tweezers experiments, the measurable quantities are the vertical extension z 
and the number of turns n imposed on the bead. Using Eq. ([3]) for n = Lk, the equation 
z = i — ip and the constitutive relation r = Mcxt/K^ where Mext is found from Eq. iQ, we 
obtain the vertical extension of the filament as a linear function of the number of turns n: 

2Ko sin^aA „ Air R 

l + —^^-)i + Xn——. 12 
ii.c! COS 2q;/ sm 2q; 



Thermal fluctuations dominantly affect the tails and make the end-to-end distance z of the 
molecule smaller than the contour length £ — £p of the tail parts, by a factor p„ic G [0, 1]: 
^ = Pwic — ^p) • This factor depends on both the pulling force Fext and the bending 
persistence length A = KQ/{kT) and can either be read off an experimental hat curve from 
the value z{n = 0) = Pwic or computed from theoretical formulas 10, l3l|]- To account for 
these thermal effects, we replace Eq. ( !T2|) with: 



Z = Pwlc 1 + 



2 Kq sin^ a 
cos 2a 



« + X Pwlc ■ „ n 
sm/a 



(13) 



One of the main features of the experimental hat curves is the linear decrease of the 
vertical extension with the number of turns. We deflne the slope q in the linear part of the 
hat curve as: 



<1 



dcf 



dn 



Pwlc— ^ • (14) 
sm 2a 

Given experimental values of Fcxt and g, Eqs. (11 01) and can be solved for R and a. Since 
q (and -Fcxt) are constant along the linear part of a hat curve, the values of R and a thus 
determined will be constant as well. As a result, the twisting moment in the molecule, given 
by Eq. ([9]), is constant, for a given experiment, along the linear reg ion of the hat curve. 



a property that has been previously reported in the literature [8|, |20| and which is a clear 
outcome of the present model. An interpretation of the fact that R and a are constant in 
the linear region of the hat curve is that each additional turn of the bead is used to convert 
a small piece of tail into plectonemes. 



Twisting moment 

The twisting moment in the molecule, which is uniform and equal to Mext at equilibrium, 
cannot be measured in magnetic tweezers experiments. However, it has been shown that 
enzyme activity such as RNA polymerase depends on the value of the twisting moment in 
DNA 12|. The value of Mext can be determined from Eq. ([9]) once R and a are known, as 
explained above. Here, we give a formula for Mext directly as function of the experimental 
slope q and the external force Fext- Indeed, using Eq. (fT4l) to eliminate R in Eqs. ([9]) and (ITOl) . 
one obtains Mcy±{q, a) and -Fext(5', a) as functions of q and a. It is then possible to ehminate 
a, which yields: 

Mext = m + (m^ + 2K, F,xt) , where m = - (15) 

47rp^ie 2q 



In the limit of small a, one can expand the functions Mcxt(q', «) and -Fcxt(q', «) prior to 
elimination of a, and this leads to a simplified formula: 

Mext ^ -^^Fext, (16) 

where, as explained above, pwic = z{n = 0)/i. As shown in Fig. HI this approximation is 
accurate when used with typical experimental values. Eq. (fT6l) provides a simple and direct 
mean of evaluating the twisting moment in magnetic tweezers experiments, based on the 
slope of the linear part of the hat curve only. Note that it should not be inferred from 
Eq. (fT6|) that Mgxt depends hnearly on Fext, as the slope q is itself a function of Fext- 

Superhelical angle limit 

It is known that the topology of contact between two impenetrable helical tubes winding 



along a common axis changes when a becomes larger than n/A (32|. The possibility of such 
a change of topology is not considered in our model (being specific to hard-core repulsion 
between tubes, it is not relevant to DNA molecules undergoing long-range electrostatic 
repulsion anyway). Nevertheless, the equilibrium solutions found here are all such that 
a < 7r/4. This upper bound has a mechanical origin, and not a geometrical one: the 
expressions for Fext in Eq. (ITOl) and for Mext in Eq. (Q both diverge at a = n/A and 
plectonemic solutions with a superhelical angle larger than n/A are unstable. 

Application to experiments 

The model is used to extract mechanical and geometrical parameters from experimental 



data. To allow comparison with previous work, we use the same data as in [2l| . These data 
are shown in Fig.O they were obtained on a 48kbp lambda phage DNA molecule in a lOmM 
phosphate buffer. 

For each curve in Fig. [2l corresponding to a given value of the external force -Fcxt, we 
extract the slope q by fitting the linear region. The superhelical variables R and a are found 
by solving Eqs. ffTOj) and f|T^ . and are plotted in Fig. [3] as a function of F^xt- The recon- 
structed values of R are in the nanometric range; they decrease with the pulling force, from 
approximately 6 to 2 times the DNA crystallographic radius in this particular experiment. 
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FIG. 2: Experimental hat curves showing the vertical extension of a lambda phage DNA 48kbp 
molecule as a function of the number of turns imposed on the magnetic bead (salt concentration 
10 mM, temperature 298 X). Experimentally measured persistence length of the molecule is ^ = 
51.35 nm. Each curve corresponds to a fixed pulling force -Fext^ 0.25, 0.33, 0.44, 0.57, 0.74, 1.10, 
1.31, 2.20, 2.95 pN. Triangles represent the fit for the slope q of the linear region. Data kindly 
provided by V. Croquette (LPS - ENS, Paris). 



At large values of the force, R is close to (and actually smaller than) the Debye length, 
3.07 nm in 10 mlvl salt, and the Ivlanning condensation radius, 3.18 nm in 10 mlvl salt 33|. 
We note that the values of R found here in the presence of a pulling force are smaller than 
(and in the same range as) in Ref. 16|] where no force is applied, which is consistent. 

The reconstructed values of the twisting moment Moxt and of the contact pressure p are 
given in Fig. HI based on the same experimental data. The values of Mext are determined 
both by Eq. ([9]) using the previously computed values of R and a, and by the approxi- 
mate formula (flGl) directly. A good agreement is obtained, which validates the proposed 
approximation. The values of Mext are also compared to those predicted by a composite 
analytical model, see Eq. (17) in Ref [2^ (this model uses effective parameters determined 
from IVIonte-Carlo simulations 1341). 
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FIG. 3: Reconstructed values of the plectonemic radius -R as a function of the puhing force, from 
the data in Fig. [2] by solving Eqs. (jlOp and (jl4p . The angle a is shown in the inset. 
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FIG. 4: Reconstructed values for the twisting moment in the molecule based on the data shown 
in Fig. [21 using the exact formula in Eq. Q (solid squares), and the small angles approximation 
in Eq. (|16p (open circles). Comparison with the prediction of the composite model in Ref. [2^ 
(curve). Contact pressure is shown in inset. 



IV. CONCLUSION 

We have shown that, under the approximation that thermal fluctuations are neglected in 
the plectonemes, one can calculate analytically the response of twisted DNA: supercoils are 
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described by a mechanically exact and self-contained model. Self-contact in the plectonemic 
region is treated with a hard-core potential; an expression for the contact pressure between 
the two dsDNA is derived. The hard-core radius is an effective parameter determined, for a 
given value of the applied force, from the slope of the linear region of the experimental curve. 
A formula for the twisting moment is proposed, as a function of the slope of the linear region 
of the experimental hat curve only. We apply this analysis to experimental data from which 
we extract the mechanical quantities: superhelical radius and angle, contact pressure and 
twisting moment. We compared these values with predictions from previous analyses, when 
available, and found that they are consistent. In future work, we shall extend the present 
model to deal with long-range interaction potentials, predict the superhelical radius, and 
utilize magnetic tweezers experiments to probe DNA-DNA electrostatic interaction. The 
present paper is a first step towards a mechanically accurate description of bare dsDNA 
subjected to tensile and torsional loads, a problem relevant to the architecture of DNA in 
the cell nucleus where proteins come into play. 

We thank V. Croquette for allowing us to use his unpublished experimental data. 
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